Geometric relation between centrality and the impact parameter in relativistic 

heavy-ion collision^! 
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We show, under general assumptions which are well satisfied in relativistic heavy-ion collisions, that 
the geometric relation of centrality c to the impact parameter b, namely c ~ wb 2 /ai ne l, holds to a 
very high accuracy for all but most peripheral collisions. More precisely, if c(N) is the centrality of 
events with the multiplicity higer than N, then b is the value of the impact parameter for which the 
average multiplicity n(b) is equal to N. The corrections to this geometric formula are of the order 
(An(b)/n(b)) 2 , where An(b) is the width of the multiplicity distribution at a given value of 6, hence 
are very small. In other words, the centrality effectively measures the impact parameter. 
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Data from relativistic heavy-ion collisions (SPS, 
RHIC) are typically categorized by introducing central- 
ity, c, defined as the percentile of events with the largest 
number of produced particles (as registered in detectors), 
or largest number of participants (as determined from 
zero-degree calorimeters) . We denote this number gener- 
ically as n. Results of measurements, such as multiplici- 
ties ]l|,p] , spectra J0-^| , the elliptic flow coefficient 
HQ, the HBT radii |gf, etc., are then presented for var- 
ious centralities. From the experimental viewpoint the 
centrality is a good, unambiguous criterion allowing to 
divide the data. On the other hand, theoreticians need to 
assign an impact parameter, b, to a given centrality. The 
impact parameter is in a sense more basic, since it de- 
termines the initial geometry of the collision and appears 
across the formalism. Theoretical calculations in heavy- 
ion physics input b in order to obtain predictions. Having 
done the calculation, the question arises as to which cen- 
trality data should the model results be compared to. For 
that purpose one typically applies the Glauber model in 
order to compute the number of wounded nucleons or 
binary collisions at a given 6, which are subsequently re- 
lated to multiplicities or number of participants P,|lO|]. 
Since these are measured in the experiment, one is able 
to identify b with c. 

In this paper we argue that such an effort is not nec- 
essary, since, under general assumptions which hold very 
well in relativistic heavy-ion collisions, we have, to a very 
high precision, the relation 



c(N) 



irb(N) 2 

Omel 



for b < R 



(1) 



where ai ne \ is the total inelastic nucleus- nucleus cross sec- 
tion, and R is of the order of the sum of the radii of the 
colliding nuclei. The centrality c(N) is the centrality of 
events with the multiplicity higher than N, while b(N) is 
the value of the impact parameter for which the average 
multiplicity n(b) is equal to N. As will be shown, Eq. 
(|l|) holds to a high accuracy for all but most peripheral 
collisions. Note that it is geometric in nature, and does 
not involve explicitly the variable n needed to categorize 



the data (multiplicities, number of participants, number 
of binary collisions, etc.). At first glance, this fact may 
seem a bit surprising. 

One can explain the geometric nature of (|l|), and the 
fact that it does not explicitly depend on n, with the 
following pedagogical example. Consider a competition 
where archers are shooting at a target of radius R, each 
of them once. The archers are very poor, such that 
they shoot randomly. They are paid accordingly to their 
aim: more central, higher reward. We are not allowed to 
watch the competition, hence do not know which spot on 
the target has been hit, but later we review the reward 
records. Suppose a large number of archers scored (here 
we take only 10 in order to write down the results ex- 
plicitly), and are ranked according to their prizes, which 
are: 100$, 100$, 50$, 50$, 50$, 10$, 10$, 10$, 10$, 10$. 
The two archers that got the highest prize (100$ in this 
case) had to hit the bull's eye. Since these are the 20% 
of all archers, and they were shooting randomly, we can 
immediately determine (neglecting the statistical error) 
the radius b of the bull's eye, since 20% is the ratio of 
the area of the bull's eye to the total area of the tar- 
get: 20% = nb 2 /(TTR 2 ). Therefore b = i?V20%. Now, 
imagine another competition is held, with all rules the 
same but the prizes differently assigned to the rings of 
the target. Suppose the ten archers got 500$, 500$, 100$, 
100$, 100$, 50$, 50$, 50$, 50$, 50$. Again, we can de- 
termine that the 20% of the highest rewards correspond 
to hitting the central spot, and can determine its radius 
exactly b as before. Note that in the determination of 
b we are not using the actual values of the rewards at 
all - the function used can be any monotonic function of 
the centrality. The rewards are only used to categorize 
the data. Once this is done, we can identify the c "most 
central" archers and determine b according to Eq. (jl]), 
irrespectively of the function used for categorizing. Our 
example can be translated into heavy-ion collisions in the 
following way: archery competition - heavy-ion experi- 
ment, archer that scored - event, rewards in competition 
I - number of participants, rewards in competition II 
multiplicity of produced particles, percentile of highest- 
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scoring archers - centrality, radii of rings on the target - 
impact parameters. 

The above example shows the essence of our argument, 
valid for the classical physics of relativistic heavy-ion col- 
lisions. There are, however, two additional features which 
need to be considered. First, a collision at a particular 
impact parameter b produces values of n which are sta- 
tistically distributed around some mean value n(b) with 
a distribution width An(b). As we will show, Eq. ([[]), 
formally valid at An(b) <C n(b), is accurate even for re- 
alistically large An(b), such as obtained from statistical 
models of particle production. Second, there are bound- 
ary effects near b ~ R — at lower values of b the inelastic 
cross section is the cross section for colliding black disks, 
whereas at the boundary the target gradually becomes 
transparent. 

We now proceed with a formal derivation. Let P(n) 
denote the probability of obtaining value n for the catego- 
rizing function (multiplicity of produced particles, num- 
ber of participants, number of binary collisions, etc.). 
For simplicity of the language we call it the multiplicity, 
bearing in mind it could be any of these quantities. The 
centrality c is defined as the cumulant of P(n), namely 



z(N) = £ P(n). 



(2) 



Thus c(N) is the probability of obtaining an event with 
multiplicity larger or equal to N. A particular value of 
multiplicity n may be collected from collisions with var- 
ious impact parameters b' , thus we can write 



c(N) =22 — 
„tVo f 



2irb'db' 



p{b')P{n\b') 



(3) 



where 2irb'db' is the area of the ring between impact pa- 
rameters b' and b' + db', the quantity p(b') is the probabil- 
ity of an event (inelastic collision) at impact parameter 
b', and P(n\b') is the conditional probability of produc- 
ing multiplicity n provided the impact parameter is b' . 
The function p(b') is unity for b' below R, and drops 
smoothly to zero at b' around R, reflecting the washed- 
out shape of the nuclear density functions at the edges. 
The interpretation of Eq. (||) is clear: the probabilities 
for hitting the ring between b 1 and b' + db', the prob- 
ability for an event to occur at b', and the probability 
to produce multiplicity n (provided the event occurred 
at b 1 ) are multiplied, as requested by the classical na- 
ture of the problem. Since we have 53 n =i P{ n W) = lj 
and, by definition, J °° 2-Kb'db' p(b') = find) we verify the 



proper normalization in Eq. (|3|), namely c(l) = 1. Fur- 
thermore, for heavy nuclei we may use the continuity 

limit > J2n=N -> In dn = So° dn 9 ( n ~ iV )- 

The function P(n\b') is not known, yet, by the statisti- 
cal nature of the particle production, and by experience 
of various models, we expect that for large values of n it 
is narrowly peaked around an average value nib'). Thus 



we begin our study by taking the limit of an infinitely- 
narrow distribution, P(n\b') = S(n — n(b')). In this case 



c(N) 



dnO(n — N) 
ii j a 

2irb'db' 



2-Kb'db' 



find 



find 

p{b')0(n(b') - N) 



p(b')5(n-n(b')) 
(4) 



Since n(b') is a monotonically decreasing function of b' , 
we have 8(n(b') - N) = 8(b(N) - b'), where b(N) is the 
solution of the equation n(b) = N. Therefore 



c(N) 



2 -^p(b>)e(b(N)-b>) 

find 

b ^27Tb'db' find^V)) 

P{b ) = : 

find find 
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where ai nc \(b(N)) is the inelastic cross section accumu- 
lated from b' < b(N). Equation (||) is a generalization 
of formula (|l|). In Ref. O] it has been quoted in the 
context of the Glauber model. We notice that although 
c and b depend implicitly on N, their relation does not 
explicitly involve N. 

We now turn to a quantitative analysis of dispersion 
effects. Assume 



P{n\b') = 



1 



: exp 



{n-n{b')f 



(6) 



An(b')V27r V 2An(b') 2 
which is a good approximation for An(b r ) < n(6'). Then 



c(N) 



2-Kb'db' 

find 



P(b') 



or 



N 



\V2An(b') 
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For small An(b') the function in curly brackets resembles 
the function 6(n(b') — N), washed out over the range 
An(b'). Thus, we introduce the function 



d(x) = \ 



erf 



(— 

\V2Ar 



9{x). 



(8) 



The integral of d(x) with a regular function f(x) can 
be expanded in even powers of An as follows (this is 
analogous in spirit to the Sommcrfcld expansion of the 
Fermi-Dirac distribution function at low temperatures): 



dx f(x)d(x) 



J=l,3,5,... 



dx' 



x=0 



with the coefficients 



a i 
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We rewrite the integral in Eq. (|7]) as J 2b' db' 
J dndb' 2 /dn, and use expansion (||) to obtain 
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c(N) 



<r iae i(b{N)) 



Cincl 



(11) 



-(An(b(N)) 2 — 



np(b(n)) db 2 {n) 
dn 



-N 



For inner b, where p(b(n)) ~ 1, the correction term is 
proportional to d 2 (b 2 (n)) / dn 2 . In the models considered 
below this quantity is proportional to 1/n 2 , and as a 
result c(N) = (?inel(b(N)) / ai ne \ + 0(An 2 /n 2 ), quantita- 
tively showing that the geometric identification (JjJ), or 
(||), is good for narrow distributions. 

In order to illustrate the above results and to ob- 
tain more detailed numerical estimates for the correc- 
tions we consider two models: a model inspired by 
the wounded-nucleon model ]l2"| , and the optical limit 
of the Glauber model |l3| for the binary collisions. A 
combination of these models has been used to explain 
the observed hadron multiplicities produced in RHIC 
f[4} |. We look at the Au + Au reaction, with the nu- 
cleus density profile pa(t) described be the standard 
Woods-Saxon function with the radius ro = (1.1274 1 / 3 — 
0.86,4~ 1/3 )fm, with A = 197, and the width parameter 
a = 0.54fm. The nucleus-nucleon thickness function is 
given by Ta(s) = dz pa(V s 2 + z 2 ), and the average 
number of wounded nucleons is 

/>00 /■2-7T 

71 (6) = 2A sds / dipT A ( \J s 2 + b 2 + 2sb cos ip) x 
Jo Jo 



1 - (1 - cjT a { S )Y 



(12) 



where, following Ref. [|14j, we take a = 40mb as 
the nucleon-nucleon inelastic cross section. The total 
nucleus-nucleus cross section obtained in this model is 
Cine] = 7.05b. The expressions for the dispersion of 
wounded nucleons produced at a given b is very com- 
plicated. Instead of computing multidimensional inte- 
grals, we explore, for our illustrative purpose, two cases: 
An ~ n, and An ~ y/H. Led by the sample numerical re- 
sults for the distributions given in Fig. 1 of Ref. |lj| , we 
take i) An = n/10, or ii) An = yffi. In Fig. 1 we show 
the results of computing c(N) according to Eqs. (|7|Jl^) 
with p(b') = 0(\/o-i n e\/ir — b'), and for the choices i) and 
ii) (dot-dashed and dashed lines, respectively). These 
are compared to nb(N) 2 / a- me \ (solid line), where b(N) is 
defined as the solution of the equation n(b) = N. The 
curves overlap within the width of the line, except for 
tiny regions at very low N, ( N < 2 ), corresponding 
to very peripheral collisions, and at large N, correspond- 
ing to b around 0. The discrepancy at large N follows 
from the fact that c(N) evaluated exactly continues to 
be non-zero till the maximum value of wounded nuclei, 
N = 2A, whereas b(N) by construction goes to zero at 
N = n(b — 0) f» 377. This effect is visible in Fig. 1 only 
for the choice i) for the widths. 

We can treat the dependence on N as parametric, and 
plot c(b(N)) vs. b(N). The results is shown in Fig. 2(a). 
Again, the model curves for c(b) for choices i) and ii) 



overlap with the curve irb 2 / Ci ne i except for very periph- 
eral (6 > 14fm) and very central (6 < 2fm) collisions. 
This behavior directly reflects the behavior of Fig. 1. 
The size of the correction of Eq. (11) is, at intermediate 
b, of the order of 10~ 3 . 




FIG. 1. Centrality in models i) and ii) (dot-dashed and 
dashed lines), and the function nb(N) 2 / CTi nc i (solid line), plot- 
ted as functions of the number of participants, N. 

As another illustrative example we consider the 
Glauber model of nucleus-nucleus collisions and analyze 
binary collisions, n = n co u. We use the optical limit of 
the model, which results in simple expressions. In this 
model 



c(N) 



2-Kb'db 1 



n—N 



A* 



n=N J ^incl 



Pg (n, b') 



(13) 



E / )[T(b>) ,]«[!- T (b')o-f - 



Cinol 



where for Pq (n, b') we have used the formula for the 
probability of the occurrence of n inelastic baryon-baryon 
collisions at an impact parameter b' O] (note that Pq 
plays the role of the product p(b')P(n\b') from the previ- 
ous discussion). Here T (b) is the nucleus-nucleus thick- 
ness function, 

/>00 />OG />OG p27T 

T(b) = ds dz A / dz B / dp x (14) 

JQ J-oc J-oo Jo 

p(\Js 2 + z 2 A )p{\J s 2 + b 2 + 2sb cos tp + z%). 



The sum in Eq. (|13|) can be carried out exactly, yielding, 
with the notation x = T (6') a, the expression 



2tt6W ( A 2 

Cinol 



(15) 



(l-x) A x N 2 F 1 (l,N-A 2 ;N+l; -) 

x — 1 



c(N) 



We perform the integration in Eq. ( |15[ ) numerically. On 
the other hand, the average number of collisions at a fixed 
value of the impact parameter b is n (6 2 ) = A 2 T (b) a. 
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Repeating the steps of the previous example results in 
the identification c(N) = c(n (b 2 )) = c(A 2 T (b) a). This 
function is compared to irb 2 / <j- mc \ in Fig. 2(b). The agree- 
ment is excellent and the two curves are indistinguishable 
except for very peripheral collisions (b > 13.5fm). 




b [fm] 

FIG. 2. (a) Centrality as a function of the impact param- 
eter for the models i) (dot-dashed line) and ii) (dashed line), 
(b) The same for the Glauber model for binary collisions 
(dashed line). Solid line shows 7rfe 2 / cr ine i . 

With the Gaussian parametrization of the thickness 
function the model can be treated analytically a bit fur- 
ther. We use T(b') = 1/(2tt/3 2 ) exp (-6' 2 /(2/? 2 )), with 
[3 — 4.6fm, which leads to a quite good approximation of 
the exact thickness function. Then 



c(N) 



27T/3 2 



Cinel 



CT /(27r/3 2 ) 



,A z -r 



E( „ )J tyy^^-yr 

2 AB 1 

^2 ~ I 'y/{^f3 2 ) ( n > 1 + A 2 - n ) > ( 16 ) 
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where I z (a, b) — B z (a, b) /B\ (a, b) , and B z (a, b) is the 
incomplete beta function. For large A 2 and small 
a/(2ir(3 2 ) the function I a /(2-np' 2 ) (n, 1 + A 2 — nj is well 
approximated by the step function 6 (A 2 a / (2ir(3 2 ) — nj. 
Replacing the sum by the integral in (RW we find the 
leading expression 



Omcl 



2irf3 2 
ABo 



N 



(17) 



On the other hand b 2 (N) = ~2(i 2 In (2n(3 2 N / {ABa)) , 
which immediately results in Eq. ([!]). Since An 2 = 
A 2 T(b)a(l - T{b)a) ~ n, the correction of Eq. (jl 



becomes -2-rt f3 / <j- mc \(l±n / N) z 



-0.2/ N, hence is very small at large N. 



-2irp 2 /a inel (l/N) 



As already mentioned, there are attempts (l4| to ex- 
plain the multiplicity of produced particles through a 
combination of the wounded nucleon model asso- 
ciated with soft processes, and production proportional 
to the number of binary nucleon- nucleon collisions, as- 
sociated with hard physics. The folding of the distribu- 
tions of wounded nucleons, n w , or number of collisions, 
"■coil, with the distribution of particles produced in an 
elementary event (by the wounded nucleon or in a sin- 
gle binary collision), may result in a broadening effect 
in the observed distribution of the multiplicity of the 
produced particles, n. However, we expect this broad- 
ening to be negligible in the ratio An/n, wich is the 
quantity controlling the accuracy of Eq. (nl). In par- 
ticular, for the wounded nucleon model m, 5l one has 



(An/n) 



2(A H ) 2 /(n w n 2 H ) + (An w /n w ) 2 , where the 



subscript h refers to the nucleon- nucleon collision. As- 
suming (Ah) 2 ~ nu, we find that the contribution from 
the first term is smaller than from the second term al- 
ready for moderately large n w , and An/n ~ An w /n w . 
This indicates that Eq. (Q) remains very accurate when 
multiplicities of produced particles are used as the cen- 
trality criterion. 
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